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1 Introduction 

In many Direct Numerical Simulations (DNS) of turbulence researchers inject 
power into the fluid at large scales and then observe how it "propagates" 
to the small scales [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. One such type of 
stirring is to take the force f(x, t) to be proportional to the projection of the 
velocity u(x, t) of the flow onto its lowest Fourier modes, while keeping the 
rate of injected external power constant. In this paper we perform a simple but 
rigorous analysis to establish bounds on the relationship between the energy 
dissipation rate (which is the same as the injected power) and the resulting 
Reynolds number. While this analysis cannot give detailed information of the 
energy spectrum, it does provide some indication of the balance of energy 
between the lower, directly forced, modes, and those excited by the cascade. 
This work is an extension of the analysis in [13, 14, 15], where the force is 
fixed (not a functional of the velocity) . 

Consider fluid in a periodic d-dimensional box of side length £. The allowed 
wave vectors k are of the form k = ^y-a, where a G 7L d is a d-dimensional 
vector with integer components. Let C be the subset of wave vectors that have 
the smallest possible wavenumber (namely, C consists of 2d elements: 
C = {±^ei, . . . , ±?j-ed\- The operator V projects the vector field 

u(x,t) = ^u(k,t)e ikx 

k 

onto the subspace spanned by the Fourier components with wave vectors in C: 
Pu(x,() = ^u(k,i)c ik ' x . (1) 

ke£ 

Obviously, V maps L 2 into L 2 vector fields; in fact, Pu is C°° in the spatial 
variables. The projection also preserves the incompressibility property. That 
is, if V • u(x, t) = 0, then V • Pu(x, t) = 0. 
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The Navier-Stokes equation is 



u+(u- V)u + -Vp = zMu + f , (2) 

p 



with f (x, t) taken in the form 



Pu(x,i) 



f(x,t) = e ■ (3) 



where || • || 2 stands for the L 2 -norm, ||"Pu(-, i)|| 2 := [/ |"Pu(x, i)| 2 d d x] 
This choice of forcing ensures that the input power is constant: 



/ 



u(x, t) • f (x, t) d d x = £ d e . (4) 



In this approach e, v and ^ are the (only) control parameters. On average, the 
power input is the viscous energy dissipation rate: 

£ :=l/u(x, i ).f(x,t)d d x = ,^M) , (5) 



where (•) stands for the long time average. The non-dimensional measure of 
energy dissipation is defined as 

V--=W^ (6) 

which is a function of Re := the Reynolds number, where U is the r.m.s. 

velocity defined by U 2 := ^"J 2 ^ , a measure of the total kinetic energy of the 
fluid. Our analysis will establish limits on the relationship between /? and Re. 

Because we will study the "low-fc" Fourier modes (i.e., modes with wave 
vectors in £), we also introduce the r.m.s. velocity V contained in these modes, 

The bounds on the dissipation (3 will be in terms of Re and the quantity 



V / "low-fc" kinetic energy of the fluid 
^ U y Total kinetic energy of the fluid 

The case p w 1 corresponds to laminar flow, when the turbulent cascade is 
inoperative. 
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2 Derivation of the bounds 



2.1 Lower bounds on the energy dissipation 

To obtain lower bounds on the energy dissipation, we proceed as usual by 
multiplying the Navier-Stokes equation (2) by u(x, t) and integrating over 
the volume of the fluid to obtain the instantaneous power balance, 

^l|uM)||^ = -H|Vu(.,*)||3+f* e , (9) 

where ||Vu(-,i)||| := / |£-, ro=1 d jUm {^ d d x. 

Now we use the facts that the lengths of wavevectors k ^ C are at least 
2ttV2/£, and that ||u(-,t) - Vu(-,t)\\ 2 2 = \\\i{-,t)\\ 2 2 - \\Vu(-,t)\\l, to derive a 
lower bound on || Vu(-, t)\\ 2 - 

||Vu(-,t)||l= / |Vu(x,i)| 2 d d x = ^^fc 2 |u(k,i)| 2 
J k 

= f ^fc 2 |u(k, i )| 2 + ^fc 2 |u(k,t)| 2 ) 

\ke£ k££ / 

^^¥fEK k ^)i 2 + 2 Ei u ( k ^)i 2 ) 

\ke£ k^£ / 



(||Pu(,t)|| 2 + 2||u(,t)-7'u(,t)|| 2 ) 



4^ 

I 2 

Att 2 

— (2||u(,t)|| 2 -||7'u(,t)|| 2 ) . (10) 



From (9) and (10) we obtain the differential inequality 

^l|u(-,t)||^<-^||u(.,t)|| 2 + £ d e , 
from which, using Gronwall's inequality, we deduce 

\ ||u(,t)|| 2 < \ ||u(,0)|| 2 e-^ + e d e^ (l - e"*^*) . (11) 



The inequality (11) implies that ||u(-, t)\\ 2 is bounded uniformly in time, which 
in turn implies that the time average of the time derivative in (9) vanishes. 
This ensures that the time-averaged power balance (assuming that the limit 
associated with the long time average exists) is indeed given by (5). 
Taking the time average of (10), we obtain the bound 
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I 2 

which in non-dimensional variables reads 

4tt 2 



-^(2-p 2 )</3. (12) 



2.2 Upper bound on the energy dissipation 

To derive an upper bound on (3, we multiply the Navier-Stokes equation (2) 
by ||p"( X t*j| 2 and integrate. The term with u gives a total time derivative, 

J \\Vu\\ 2 \\Vu\\ 2 J dt y ' 2 dt 

For the viscosity term we obtain (integrating by parts) 



v J (An) 



Vu , llVPull 2 4tt 2 



|Pu|| 2 WVuh ' t 2 



2 



while the forcing term gives £ d e/\\Vu(-, t)|| 2 (cf. (4)). 

To estimate the inertial term, we introduce temporarily the notation 
p(x, t) := Vu(x, t). We will make use of the uniform (in x and t) estimate 

\d jPm (x,t)\ ^ \kj\ \u m (Kt)\ _2*_ E ke£ |fl(M)l 2i:Vd 
\\p(;t)h ~tc H p (-'*)ll 2 " ^-ecH*,t)V ~ e + i ' [ ' 

Then the inertial term may be bounded (we use V • p = 0): 
f [(u • V)u] • d d x = [u ■ • ud d x 

J \\P\\2 J || E> || 2 

27T\/d /" 2 d 27T\/d 2 

This estimate, however, is obviously not going to be tight for small Re, when 
the flow is not very turbulent. To improve this estimate so that it take into 
account the fact that for small Re the energy does not "propagate" much 
from the large to the small wavenumbers, we split the velocity u into a "low- 
k" component, Vu, and a "high-fc" one, u-Pu. We will still use the uniform 
estimate (13) as well as the inequality ab < ^(za 2 + jb 2 ) which holds for any 
z > 0: 
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/ [(u . V )u] • ^ d*x| . |/[p + ,u- p)] • • [p + (u - p)] d'x 



£l + (d/2) 

2ny/d 



. -^—^ J (2|p||u-p| + |u-p| 2 ) ,r'x 

< J^r/Npi 2 +a+i) i" -pi 2 ] ^ 

< 1(1 + 1) Ml +(z-^-l)\\P\\l] . (15) 
Putting together (14) and (15), we find 

fd e 1 < l d||^u(-,t)|| 2 

' \\Vu{;t)h ~ 2 dt 

+ $m ^{H\l (i + 1) Hi 2 . + (* - \ i) IIp»H3} 

+ ^||Pu(,i)|| 2 . (16) 

Now take the time average of all terms in the above inequality. First note 
that the average of the time derivative of ||'Pu(-, t)|| 2 gives zero due to the 
boundedness of ||7^'u(-, i)|| 2 (which follows from the boundcdncss of ||u(-,i)|| 2 ; 
sec (11)). To estimate the other terms, we use Jensen's inequality: if a function 
9 is convex and (•) stands for averaging, then (9 o g) > 9 ((g)) for any real- 
valued function g. Applying this inequality to the case g(t) = ||"Pu(-, i)|| 2 and 
the convex function 9(t) = t 2 , we obtain (same as Cauchy-Schwarz) 

(\\Vu\\ 2 )<^(\\ru\\ 2 )=£ d / 2 V . 

On the other hand, if we take 9{t) — j for t > 0, we deduce 

1 \ 1 1 1 

> _ __ > 



\Vu\\ 2 / ~ (\\Puh) ~ y/{\\Pu\\%) f d/2 V 
Plugging these estimates into (16), we obtain 

e < + ^ min {U*V, (± + l) U 2 V + (z - j - l)V 3 } . 

In terms of the non-dimensional energy dissipation rate (6), we can rewrite 
this inequality in the form 

4_2 

(3< — p 2 + 27rVd0(p,z) , (17) 

where we have introduced the function 

<Xp, z) := min {p, (I + l) p + [z - \ - l) p 3 } . (18) 
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2.3 Compatibility of the lower and upper bounds on /3 

Assembling the lower and upper bounds (12) and (17), we have 

^(2-p 2 )</3<^p 2 + 2W^(p,z) • (19) 

The compatibility of the two bounds on f3 imposes restrictions on the 
allowed range of p, namely, p should satisfy the inequality 

P 2 + ^%(P^)-1>0. (20) 

In the interval p <E [0,1], this inequality is satisfied for p e [p m ; n (Re, z), 1], 
where p m ; n (Re, z) w J^ Re for large Re. Clearly, the lower bound on the range 
of p is more meaningful for smaller Re. 



2.4 Optimizing the upper bound on /3 



Since we do not have a priori control over p, we will derive an upper bound 
for P by maximizing the upper bound in (19) over p, after which we use the 
freedom in the choice of the parameter z > to minimize for any given Re, 
which results in 



(3 < min max 

z>0 pe[p min (Re,z),l] 



2 



^p 2 + 2Wd0(p,z) 



(21) 



Although this procedure is not difficult to implement numerically we will 
analyze only the case of high Re where the analysis can be carried out an- 
alytically. First notice that for high Re, the lower bound p m ; n (Re, z) is very 
small, so the maximization over p can be taken in the entire interval [0, 1]. 
Thus <f>{p, z) < <j,*{z) := max pe[0!l] <f>(p,z) = (l + z-z 2 )" 1 / 2 for z G [0, 1± ^)- 
Since for large Re the Re-independent term in the right-hand side of (21) is 
dominating, we have the high- Re estimate 



(3 < min 



2 



^<t>*(z) 2 + 2nVd^(z) 



16tt 2 4wVd 



+ —7^ ■ (22) 



5Re V5 



At high Re, the value or p maximizing <fi(p, z) is -?=. We remark that it is not 
difficult to prove that the upper bound (22) is optimal (i.e., coincides with 
( 21 )) forRe ^57b- 

3 Discussion 

In dimension 3, the scaling of the upper bound is in accord with conventional 
turbulence theory: at high Re, e ~ ^4- is independent of the molecular viscos- 
ity. For the type of forcing considered here, we find /3 < Air J | w 9.7339 . . .. 
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Fig. 1. Upper and lower bounds on [3 (solid and dashed lines, resp.) 



A plot of the bounds is presented in Figure 1. At low Re, the upper and lower 
bounds converge to each other. While it is difficult to compare these bounds 
quantitatively with DNS results, we note from [7] that at high Re, values 
of (3 are typically around 1. Hence, our rigorous analysis, while yielding the 
expected scaling, overestimates the constants by about an order of magnitude. 

In the 3-dimensional case, if we assume that the cascade is Kolmogorov, 
i.e., the spectral density of the energy is given by E^{k) — Ce 2 / 3 fc~ 5 / 3 , we 
can estimate the "Kolmogorov" value pK as follows: 

f°° 27T /27r\ 

Ekin, total ~ J E K (k)dk , £ki ni low fc ~ — E K ^ — J , 



which gives pK 



Plugging this value in (21) and minimizing over z, we 



obtain the (approximate) estimate 



26.3 
"Re~ 



+ 8.9 



which gives a slight improvement compared with the bounds (22). 
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